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We introduce axiomatically a complete thermodynamic formalism for a single macromolecule, 
either with or without detailed balance, in an isothermal ambient fluid based on its stochastic 
dynamics. With detailed balance, the novel theory yields mesoscopic, nonequilibrium generalizations 
for entropy (T t ) and free energy (^t) of the macromolecule. T t and "ft fluctuate. Expectation 
(d/dt)E [ift] < 0, "=" holds if and only if the macromolecule is at thermal equilibrium, in which we 
show Tt still fluctuates but iPt is a constant. The entropy fluctuation a la L.D. Landau, E [(AT t ) 2 ] , 
precisely matches the fluctuation in the internal energy, which in turn equals the fluctuation in heat 
dissipation. As a generalization of Clausius' classic result, the dynamic fluctuations in the entropy 
and energy of the macromolecule are exactly compensated at thermal equilibrium. For systems with 
detailed balance, Helmholtz free energy is shown to be the potential of Onsager's thermodynamic 
force. 
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In recent years, the stochastic model for overdamped Newtonian Brownian dynamics in a force field, 

EdX t = F(X t )dt + TdB f , (1) 

has found an increasing number of applications to a host of macromolecular processes in equilibrium and more 
importantly nonequilibrium steady-state (NESS) [0J^]. These models are generalizations of the classical theory of 
polymer dynamics in which F usually is conservative || . In addition to obtaining the stochastic dynamics from the 
model, however, it becomes clear to us that an axiomatic isothermal thermodynamic formalism can be developed 
based solely on the stochastic differential equation (jl|) in which X t represents the coordinates of the "atoms" in the 
macromolecule, S is a constant, positive definite hydrodynamic interaction matrix, F is the force not necessarily 
conservative [@, and TdB t is a white noise representing the collisions between the macromolecule and the solvent: 
rr T = 2fceTS according to Einstein's relation and defines the temperature of the isothermal system. The solution 
to (|l|), X t , is stochastic, whose probability density function satisfies the Fokker-Planck equation 

^ = V- (^AVP-E- 1 F(x)Py (A = S-!rr T S- T = 2k B T~- 1 ), (2) 

where P(x,t) is the probability of the macromolecule being in conformation x at time t: X t = x. In this paper, we 
show how a complete, mesoscopic thermodynamic theory can be formulated based on (||), which we call stochastic 
macromolecular mechanics (SM3) |^]. We then apply this novel formalism to further elucidate a classic observation 
in the equilibrium thermodynamics of proteins: the dynamic origin of entropy-enthalpy compensation [^J. Naturally, 
the enthalpy under constant pressure is equivalent to the internal energy in our formalism. The significance of 
the nonequilibrium steady-state obtained from this formalism || and its applications to free energy transduction in 
biological macromolecules, e.g., a protein molecule converting chemical potential into mechanical work, have been 
discussed elsewhere ■ 

It is important to point out that even though there is already a vast literature on stochastic models based on 
Fokker-Planck equation (Q) , it is not known that this approach rigorously encompasses a comprehensive statistical 
thermodynamics. Furthermore, it is not known whether a thermodynamics exists for individual macromolecules. In 
this paper, we demonstrate the logical relation between models based on equations (Q) and (|^) and the principles 
of equilibrium and nonequilibrium thermodynamics. A complete statistical thermodynamics for isothermal systems 
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in equilibrium, nonequilibrium transient and NESS, as an complementary and alternative to Boltzmann's approach, 
emerges. 

Stochastic Macromolecular Mechanics. Following Lebowitz and Spohn (|||, we first introduce the instanta- 
neous heat dissipation Wt'. 

dW t = F(X t ) o dX t = F(X t ) ■ dX t + ^dX t VFdX t (3) 

where o and • denote integrations in the Stratonovich and Ito sense, respectively |^|. According to this definition, the 
heat dissipation, dW t , is equal to the work done by the system: the product of force F and displacement dX. This is 
the law of energy conservation. Using the expression (|l|) it is easy to show that the mean rate of the heat dissipation 
(HDR, h d ) ffl 



d f 
h d = J t E[W t ] = / F(x)-J(x, 



t)dx (4) 



where 



J{x, t) = -^AVP{x, t) + E~ l F{x)P{x, t) (5) 

is the probability flux in (|J): dP(x,t)/dt = —V • J. In mechanical terms, the mean heat dissipation rate is the 
product of the molecular force and the flux. 

The dynamic of Wt is itself Brownian- motion like with hd as its mean rate || . Its fluctuations can be characterized 
by a "heat diffusion coefficient" : 



E 



{dW t f 



(II 



= a ,r/ W E- W (,^ (6) 



which has a dimension of [energy 2 ] / [time] . The statistical properties of Wt have been extensively explored in connection 
to the Gallavotti-Cohen symmetry and fluctuation theorem |h|||,[|. 

The Onsager's thermodynamic force is different from the mechanical force F{x) |nj|. In terms of (||), we introduce 
a second thermodynamic quantity, n(a;,t), the thermodynamic force 12 



n(x,t)= F(x)--EAV In P(x,t), (7) 

and a third thermodynamic quantity: the entropy S according to the well-known formula: 

S = -k B J P(x, t) In P(x, t)dx (8) 

where P{x 7 1) is the solution to In terms of (0) and @ we have the increase of the entropy at constant temperature 
T 



TS = kgT J (InP + 1)V • J dx 

= - J (fcsTVlnP -F)-Jdx-J F J 



dx 



II • Jdx - h d (9) 

= e p - hd- 

in which we identify, following Onsager, J II • Jdx with the entropy production rate (EPR, e p ) |p"3f . It can also be 
rewritten as 



J QaV1hP(x,<) -E- 1 F(x)\ E^A^]nP(x,t)-E- 1 F(x)\p{x,t)dx. 



(10) 



which is always non-negative. This is the second law of thermodynamics. (Q) is valid for all isothermal Brownian 
dynamical systems, with or without detailed balance, in stationary state and in transit process. It encompasses both 
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the first and second laws of thermodynamics. In fact, it makes the second law quantitative by providing a rate for 
entropy increase. In a time independent stationary state, the S = in (^|), and the entropy production e p is balanced 
by the heat dissipation hd- This is the general case for an isothermal NESS. 

Eq. [Io| also indicates that e p equals zero if and only if F = (SA/2)VlnP = V In P/k B T. That is the force 
field F has to be conservative with an internal potential energy: F = — VJ7. For system with the potential, also 
known as detailed balance JTJ], the stationary solution to tffl) is P = Z~ x e~ u l kuT where the normalization factor 
Z is the partition function in Gibbsian equilibrium statistical mechanics (isothermal canonical ensemble). It can be 
mathematically shown that e p = is a sufficient and necessary condition for the stationary stochastic process X t 
defined by ([!]) being time reversible (|T5| . Therefore, time reversibility, detailed balance, and zero entropy production 
are equivalent with an equilibrium, fllpf 

For systems satisfying the potential condition F = — VJ7, the thermodynamic force II also has a potential: II = 
-V* where ^(x) = U(x) + k B T\nP. We note that the expectation of ^ t = ^(X t ): 



E [* t ] = E [U(x) + k B T\nP{x)] = J P(x, t)U(x)dx - TS 



(11) 



which in fact is the Helmholtz free energy ! The first term in ( |TT|) is the mean internal energy. More importantly, 
it is easy to show that 



E 



rS-=-e p <0. (12) 



In an isothermal system, the Helmholtz free energy decreases and reaches its minimum at the equilibrium: —k B T\nZ. 
This is precisely the statement of second law of thermodynamics for an isothermal system. 
Finally, with the potential condition, we have 

- dU(X t ) = -VU(Xt) ■ dX t - ^dX t ■ VVU(X t ) ■ dX t 

= F(X t )-dX t + ^dXfVF(X t )-dX t . (13) 

Comparing (13) with (J3J) , we see that the heat dissipation dW t = —dU(X t ), the internal energy fluctuation. Hence, 
Wt — —U(Xt) is stationary and its expectation and variance are the internal energy and heat capacity (C v ) of a 
single macromolecule at thermal equilibrium. 



The Matching Entropy and Energy Fluctuations. We now focus on systems with detailed balance. The 
above thermodynamic formalism suggests the fluctuating U t = U(X t ) as a mesoscopic, nonequilibrium generalization 
of internal energy of a macromolecule in an isothermal aqueous solution. Then its expectation E[Ut] — J U (x)P(x)dx 
which equals to the standard internal energy in thermodynamics. Similarly in the same spirit, we can introduce 
fluctuating entropy T t = — kg In P(X t ), which can be viewed as the mesoscopic, nonequilibrium generalization of 
entropy. It is important to point out that this definition is consistent with the Boltzmann's microscopic entropy based 
on the volume of the phase space. In our case P(X t ) is the probabilistic measure of the phase space. The expectation 
E[T t ] is the Gibbs entropy in Eq. ||. 

We now show that the mesoscopic generalizations immediately lead to an interesting thermodynamic result in 
equilibrium. We note that while Ut and T t are fluctuating in equilibrium, their difference, 

U t - TT t = U(Xt) + k B T\nP{X t ) = -k B T\nZ (14) 

the equilibrium free energy, however, is not fluctuating. We can further compute the fluctuations in the mesoscopic 
entropy and internal energy: 

T 2 £[(AT t ) 2 ] =E[(AU t ) 2 ] 

k B T 2 -^E{U t } = k B T 2 C v . (15) 



Furthermore, 



TE [AT t AU t ] — E [(AC/ t ) 2 ] . (16) 
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Therefore, the fluctuations of the mesoscopic internal energy Ut and entropy T t are perfectly correlated. They 
compensate in the dynamical fluctuationsi of a macromolecule. 



Entropy Fluctuation and Its Mesoscopic Interpretations. The relation between entropy fluctuation and 
heat capacity in Eq. [l5| was known to L.D. Landau who also advocated the concept of entropy fluctuation |T^]. ft 
has been a difficulty concept to many who consider entropy to be a functional of the distribution [[l9| . Here we offer 
a more plausible interpretation of the concept based on our mesoscopic view from the previous section. 

Let's consider ./V identical, independent macromolecules in the aqueous solution, (Xi, X2, X n), each with its 
own stochastic dynamic equation (|l|). The concentration of the number of molecules in conformation x can then be 
defined as 

N 

C x , t = J2^~ x j)- (17) 

3=1 

Note that since the X's are stochastic, the concentration C x .t fluctuates. However, the classic statistical mechanics is 
only concerned with the most probable C x ,t since the relative fluctuation in C x ,t is insignificant in the thermodynamic 
limit when N — ^ 00. Nevertheless, the C X: t fluctuates. 

The macroscopic entropy is defined as a functional of the density function C Xt f Therefore, it fluctuates with the 
C x .t- We now show that this fluctuation is indeed the mesoscopic fluctuation introduced in the formalism for single 
macromolecules. We consider a general thermodynamic quantity 

Q{x)C x . t dx. (18) 

If one had neglected the fluctuation in C x j, a fluctuation in q would be inconceivable. 
Noting Eq. [IT], the expectation of q is readily computed: 

N 



E[q}= f Q{x)dx^E[8{x- Xj)] 
3=1 

= N J Q{x)P{x)dx = NE [Q{X t )\ 



where the joint probability for X\, X2, Xn, P(x%, X2, ...xpf) = P(xi)P(x2)---P(xn) since the macromolecules 
are assumed to be independent in the solution. Similarly, the variance in q: 



N r 

VAR[q] = VAR / Q(x)S(x - Xj)d; 



3=1 

N 

= J2VAR[Q(X 1 )} 

i=i 

= NVAR [Q{X t )\ = NVAR [Q t ] 

in which Q t = Q(X t ) pc[ |. Therefore, the fluctuations in q due to fluctuating in the distribution function C Xi t, 
is exactly the mesoscopic fluctuation we have introduced in the previous section for single macromolecules. The 
mesoscopic view, however, clearly indicates a stochastic dynamic origin of these fluctuations. 



Entropy and Heat. It is a classic result of Clausius that entropy change equals to heat dissipation in isother- 
mal quasi-static processes. Our result suggests that both concepts can be generalized dynamically to isothermal 
equilibrium, and their fluctuations are indeed equal, as we have shown. In our formalism, the energy conservation 
is instantaneous, hence dWt = —dUt- The free energy, on the other hand, has to be constant over the entire con- 
formational space in an equilibrium, ^>(x) = — fc^TlnZ, while T(a;) and U(x) are not. This is a demonstration 
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for the concept of local equilibrium, which is essential in the general theory on nonequilibrium thermodynamics pl| . 
Therefore, TdT t — dUt — —dWt] and dUt/dT t = T, all instantaneously. 



The Rate of e p . The e p defined in (jlfj) is instantaneous and time-dependent in a nonequilibrium transient. Thus, 
one can further compute its time derivative (de p /dt), the change in entropy production rate |Q. With some algebra 
we have |22| 

^ = -2k B T J (v • jf p- 1 (x)dx - J (n T s- x n) (V • J) dx. (19) 

Near a NESS, V • J is small and the second term has the leading order. However, with detailed balance and near an 
equilibrium, II is also small; hence the first term, which is negative, becomes the leading term. Therefore, near an 
equilibrium, the e^monotonically approaches to zero. This is the Glansdorff-Prigogine's principle of minimal entropy 



production rate [Elf] . Near a NESS, however, the second term dominates, and (19) is not necessarily negative 



Discussions and Summary. While our thermodynamic formalism is strictly for isothermal processes, it can be 
used to compare a system under different temperature. In a thermal equilibrium, as we have shown, the fluctuation 
in internal energy Ut = U(X t ) is directly related to the heat capacity: C v — dE[Ut]/dT. Such relations, if any, in 
nonequilibrium steady-state have not been explored. Furthermore, one can easily introduce an external force into the 
Eq. |l| to represent the pressure, and thus we expect a thermodynamic formalism for isobaric systems can be developed 
in parallel. 

It is also noted that there is a difference between the mesoscopic energy U t = U(X t ) and entropy T t = 
— k B In P{X t ). While the former can be computed along a stochastic trajectory, the latter can not until the probability 
distribution function P(x, t) is known. This difference reflects the fundamental difference between the two physical 
quantities, energy and entropy. The former is local while the latter is non-local due to circular balance: one can not 
know the value of entropy of a state until knowing how likely it occurrs in comparison with othe states. It is possible 
to formally express the P{x.t) in temrs of a path integral p3f. Maes also suggested a space-time approach to the 
problem [Q. In a thermodynamic equilibrium, however, P{x)ca.n be determined locally up to a normalization factor 
due to detailed balance. Hence, dT t = TdUt- 

In summary, 

EdX t = F(X t )dt + TdB tl (conformational dynamics) 

dW t = F(X t ) o dX t , (heat dissipation) 

Tf = -k B \nP(X t ,t), (entropy) 
U t = F(X t ) + TVT(X t ,t) (thermodynamic force) 

is a complete set of equations which provides the stochastic dynamics of a macromolecule, its heat dissipation, its 
entropy (Boltzmann), and its thermodynamic driving force (Onsager). With this set of equations, one can compute 
hd = (d/dt)E[Wt] and e p (> 0) from entropy balance T(d/dt)E[Y t ] — e p — hd- If a system is detail balanced, i.e., F 
= -VU. Then W t = U t = U(X t ), and II = -V*(X t ) where *(X t ) = % = U t - TT t is free energy: {d/dt)E[^ t ] 
=-e p < and reaches its minimum — kgT hi Z at equilibrium. In the equilibrium, the probability distribution for X t 
is z- 1 e- u ( x V kaT , U t - TT t = ~k B T\nZ, and £[(AT t ) 2 ] = (l/T 2 )E[(AU t ) 2 } = k B C v , the heat capacity. 
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